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Abstract — Recently heterogeneous base station structure has 
been adopted in cellular systems to enhance system throughput 
and coverage. In this paper, the uplink coverage probability for 
the heterogeneous cellular systems is analyzed and derived in 
closed-form. The randomness on the locations and number of 
mobile users is taken into account in the analysis. Based on the 
analytical results, the impacts of various system parameters on 
the uplink performance are investigated in detail. The correctness 
of the analytical results is also verified by simulation results. 
These analytical results can thus serve as a guidance for system 
design without the need of time consuming simulations. 

I. Introduction 

The ever increasing demands for high data rate and low 
outage probability have driven researchers around the world 
to keep developing new technologies for wireless communi- 
cations. Heterogeneous cellular system with various tiers of 
base stations (BSs) is one of those new technologies proposed 
recently to increase coverage and system throughput A 
practical heterogenous cellular system as shown in Fig. 1 
involves two tiers of BSs: macro BSs and pico BSs. With 
the deployment of pico BSs, the system throughput in hot- 
spots and/or the coverage in black holes can be significantly 
improved. Usually pico BSs transmit at a lower power compar- 
ing to macro BSs so that significant interference to the existing 
macro BSs is avoided. However, this low power transmission 
of pico BSs limits mobile users from connecting to them. In 
order to compensate this low power transmission, a power 
offset is generally adopted to favor the selection of pico 
base stations |[2|. This scheme is called cell range expansion. 
Although the benefit of heterogeneous cellular systems with 
range expansion has been justified by some simulations, a 
theoretical analysis for such systems which can provide a 
guidance to system design is still lacking. 

In this paper, the uplink coverage probability for a hetero- 
geneous cellular system with two tiers of BSs and range ex- 
pansion is analyzed. Due to high mobility of mobile users and 
random deployment of pico BSs, the locations and numbers 
of pico BSs and mobile users are essentially random. Such 
randomness is taken into account in the performance analysis. 
This is significantly different from most of the analysis for 
cellular systems which only consider fixed topology of BSs 
and mobile users O. Based on the theory of stochastic 
geometry which is a powerful tool in modeling the randomness 
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Fig. 1. A heterogenous cellular system with two 
tiers of base staions and range expansion. 



of nodes and has been adopted recently to analyze a number of 
wireless systems in Bl- lfTol . the uplink coverage probability 
is derived in closed-form. The correctness of the analytical 
results is verified by simulation results and the impacts of 
various system parameters on the coverage probability are 
discussed in detail. Specifically, the benefits of deployment of 
pico BSs and range expansion are theoretically justified by the 
analytical results. This analytical result provides an efficient 
way to evaluate the performance under various settings and 
thus provides a guidance for system design without the need 
of time consuming simulations. 

The rest of this paper is organized as follows. System model 
is introduced in Section 11. The probability distributions of 
the distances between a mobile user and its closest macro 
and pico BSs are analyzed in Section 111. With the probability 
distributions of the distances, the uplink coverage probability 
is then derived in closed-form in Section IV. In Section V, 
the analytical results are verified by simulation results and 
the impacts of various system parameters including power 
offset and intensity of pico BSs are discussed in detail. Finally, 
conclusions are drawn in Section VI. 

II. System Model 

A heterogeneous cellular system with two tiers of base 
stations, i.e., macro and pico base stations, is considered in this 
paper. Generally macro BSs are placed in a planned manner 
and distributed deterministically. A widely used hexagonal 
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Fig. 2. A hexagonal cell with two tiers of BSs 



model is thus adopted here to model the cell with macro 
BSs distributed at the center and vertexes as shown in Fig. 
2. The radius of the hexagonal cell is denoted as d and the 
transmission power of macro BSs is Pi. On the other hand, 
the pico BSs are usually employed to increase the system 
throughput in hot-spots or to improve the coverage in black 
holes. The number of pico BSs is variable and their locations 
are also irregular It is thus reasonable to model them as a 
stationary Poisson point process $2 with intensity of A2. The 
pico BSs generally transmit at a low transmission power P2 
with P2 < Pi- Similarly, due to high mobility of mobile users, 
the users are assumed to form another stationary Poisson point 
process with intensity of Aq. Without loss of generality, the 
mobile users are assumed to transmit with the same power as 
Po- 

For an arbitrary transmitter located at z, the channel between 
the transmitter and its destination located at x considers both 
the small scale Rayleigh fading effect and path loss effect. The 
channel power gain is then modeled as hxzl{\x — z\) where 
hxz captures the Rayleigh fading effect and is modeled as an 
exponential random variable with unit mean, while — z\) 
denotes the path loss and is given by l{\x — z\) ~ \\x — 2|| " 
with path loss exponent of a. Since the path loss exponent 
is larger than 2 in most cases ifTTl . a > 2 is adopted here. 
The signal power received at the destination can then be 
formulated as hxzPi/\\x — zW"' , i = 0,1,2. Consequently, 
the uplink signal-to-interference-and-noise ratio (SINR) for a 
mobile located at z and its associated BS located at x is given 

by 

hxzPo\\x - z\\ 



SINR^ = 



J2 ^xz 

z,e<s>°/{z} 
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where ^ hxziPo\\x — Zi\\ " represents the interference 

z,e<s>o/{z} 

from other mobile users and cr^ is the noise power. It is 



generally assumed that the signal can be successfully de- 
tected/transmitted when the SINR is larger than a predermined 
threshold, which is a system parameter depending on modu- 
lation and coding schemes. 

In the considered cellular system, each mobile can be 
served by either a macro or a pico BS. In practice, BS 
selection is performed based on the strength of the received 
downlink reference signals from the base stations. In order to 
avoid frequent handover due to temporary fading effect, the 
downlink reference signal strength averaged over a certain time 
period is taken as a metric for BS selection. It means only 
path loss effect will be taken into account in BS selection 
and the user will connect to its closest base station if BSs 
transmit at the same power Denoting distances between a 
typical user and its closest macro and pico BSs as ri and 
r2 respectively, the typical user will connect to the closest 
macro BS if Piri~" > P2''2 otherwise it will connect to 
the closest pico BS. In order to alleviate the load of macro 
BSs and expand the range of pico BSs, a practical range 
expansion method jl] is adopted here and the BS selection 
is conducted based on an offset power criterion. Specifically, 
a typical user will connect to its closest macro BS only if 
Piri~" > P2Ar2^", where A is the power offset, otherwise 
it will connect to its closest pico BS. Since the distances of ri 
and r2 are key factors for BS selection and uplink performance 
analysis, their distributions will be derived first in the next 
section. 

III. THE DISTANCES TO THE CLOSEST BASE STATIONS 

As Stated above, the macro BSs are located at the center 
and vertexes of a regular hexagon as shown in Fig. 2. Due 
to the symmetry of the hexagon, only an equilateral triangle 
with macro BSs located at its vertexes (as shown in Fig. 3) is 
required to consider for the analysis of the distance ri between 
a mobile user and its closest macro BS. Since mobile users 
form a stationary Poisson point process, they are uniformly 
distributed in the considered area of the triangle. Therefore, 
the distance ri between a mobile user and its closest macro BS 
can not be larger than dj \/i and the cumulative distribution 
function (CDF) of ri equals 



Fi(r)=P(ri <r) = |, 



(2) 



where S is the area of the triangle OED in Fig. 3 and A is 
the area of the intersection of the triangle OED and the union 
of three circles centered at the vertices of OED with radius r. 

Clearly, the area of the triangle OED (i.e., S) is \fid^ j ^ 
and the sectors OLQ, EKN and DPM will not intersect with 
each other when their radiuses are less than rf/2. It is easy to 
obtain the area of the union of the three sectors (i.e.. A) as 
A = 3 X ^ when r < d/2. It follows that in this case, the 
probability that ri < r is given by 



P (n < r) 
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Fig. 3. Distribution analysis of ri. 



When d/2 < r < the derivation of the intersecting 

area of A is a little bit more complicated. Considering the 
overlapping area B of the three sectors OLQ, EKN and DPM 
as the gray area shown in Fig. 3, after some tedious geometric 
calculations, we have the overlapping area as 



B = 3x 
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It follows that when d/2 < r < d/^/i, the area of the union 
of the three sectors (i.e.. A) is given by 



A = 7rr72 - 3 



12 d d I 
X I r arccos ^ 2 



(5) 



Therefore, the probability of ri < r when d/2 < r < d/\/3 
can be expressed as 

2TT /r \ 2 



P (ri < r) 
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(6) 



Finally, the CDF of ri can be written as (|7| as shown on the 
top of next page. 

With respect to the distance of r2 between a mobile user 
and its closest pico BS, the probability that r-2 < r equals the 
probability that there are no pico BSs in the circle centered 
at the user with radius r. The CDF of r2 can be derived as 
follows 

F2(r) = P(r2 < r) = l-P{$2(B(z, r)) = 0} = l-e^^^'^''', 

(8) 

where z is the location of the mobile user, B{z,r) is a circle 
centered at z with a radius of r and $^(i?(z,r)) denotes the 
number of pico BSs in B{z, r). 

In the following section, these distributions will be used to 
analyze the uplink performance of the considered heteroge- 
nous cellular systems with macro and pico BSs. 

IV. UPLINK PERFORMANCE ANALYSIS 

In this section, the coverage probability of the considered 
uplink heterogeneous cellular system, which is defined as the 
probability that a mobile user can successfully communicate 
with its associated BS, is analyzed. In fact, this coverage 



probability is the complementary of the outage probability of 
a mobile user and has been widely adopted as a performance 
metric for system evaluation ID, ||6l, ITSl . Since each mobile 
user could access to either its closest macro BS or to its closest 
pico BS, the uplink coverage probability Pc can be written as 
the sum of two probabilities: the probability Pmsuc that the 
user connects to its closest macro BS and conducts a successful 
transmission and the probability Ppsuc that a user connects to 
its closest pico BS and conducts a successful transmission, 
namely 



Pc 



(9) 



Here a successful transmission means the signal can be suc- 
cessfully decoded by the receiver, i.e., the corresponding SINR 
is larger than a predefined threshold. In the following, the two 
probabilities P^ 

According to the definition, P^ 



msuc and Ppsuc will be derived individually. 

can be formulated as 



P — E 

msuc — ^ri.r2 



1 r2 > 



P2A 
Pi 



ri 1 (71 > Ti) 



(10) 



where Ti is the threshold above which macro BSs could 
decode successfully, and 71 is the received SINR at the user's 
closest macro BS and equals 



71 



hiPaljri) 
h+N 



(11) 



where N is the additive noise power, hi denotes the exponen- 
tially distributed channel power gain between the user and its 
closest macro BS, and Ii represents the received interference 
as 

/l= h,z^Poli\x-Z^\)■ (12) 

In ( fT2] ). z represents the location of the typical user, 4>°/ {z} is 
the set formed by the points of except 2, and x denotes the 
location of the closest macro BS concerned. Since the event 

^'2 > (^P^) " independent with the event 71 > Ti, the 
probability corresponding to the closest macro BS in ( fTOl i can 
be rewritten as 



P 7-2 > 



P2A 



ri\ri P(7i > Ti|ri) 



Based on (fTTI ). P (71 > Ti\ri) can be reformulated as 



3) 



P(7i>ri|ri)=PUi> 



Ti ih + N) , 

Pol in) 



Et 



ri(Ji + JV) 

Pol(ri) 



(14) 



Since cellular networks are generally interference dominated, 
the additive noise N can be neglected. Based on (fT2b and 
Slivnyak Theorem regarding to the interference lfT2l . we can 

get 



P(7i >ri|ri) = e- 



^in(2,r/Q) 



(15) 



Meanwhile, we have 



(16) 



Putting ( fTsT i and ( fT6b into ( flj] ), we have 







P9 A \ — \ 2 



(17) 

Using integration by parts and with the CDF of ri in (|7]i, (fTTI i 
can be derived as 







n/3 



— - — ri'^arccos e ^dri, (18) 



where if equals 



a sm (27r/Q!) \ Pi J 

Based on 2.322 in ITSl . the second term in the right hand side 
of (fTsT l can be rewritten as 



AKtt 2 -/fr^ , 27r / -K!^\ 27r _^<i^ 
■rfe irfri^ I 1 - e I " a 



73^2 A' 



3V3 



(20) 

Similarly, with 3.382 in lfT3l . the third term in the right hand 
side of ( fTsT i can be derived as 



j ^^12r2-3rf2rie-^'-?dri 



ifd2 I U' / I 2' 12 



(21) 



where r(.,.) is the upper incomplete gamma function de- 

oo 

fined as r(c, x) — J e^^f^^^dt. With the properties that 



T{c+l,x) = cr (c, x) + x'^e-^ andr(i,x2) ^ y/^ - 
spa (l — 2(5 (v^a;)), where Q(.) is the Q function, the third 
term can be further given as 



v/3 

y ^y^l2r2-3d2rie-^'^?dri 
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(22) 



Substituting (|20| | and (|22] | into dTSI l and after tedious 
computation, we finally can get the probability cor- 
responding to the closest macro BS, Pmsuc, as 
as shown on the top of next page, where P'g 



Va — I X2iT[ 

J rf arccos -^e 



dri. 



Similarly, the probability that a user connects to its closest 
pico BS and conducts a successful transmission is given by 



P — F; 

psuc — -L-'ri ,r2 



1 ^2 < 



Pi 



ri 1 (72 > T2) 



(24) 

where T2 is the threshold above which the pico BS can 
successfully decode the user's signal, and 72 is the received 
SINR at the user's closest pico BS and defined as 



72 



h2Pol{r2) 



(25) 



h+N 

In dZSl ), /i2 denotes the channel power gain between the user 
and its closest pico BS, and I2 is the interference seen at the 
pico BS and defined as 



^2 = ^ h^:z,Pol {\X - Zi\), 



(26) 



where z represents the location of the typical user and (l>°/{z} 
is the set formed by the points of except z, while x denotes 
the location of the closest pico BS. Under the interference 
dominated cellular system and based on Slivnyak Theorem 
1121 . we can get 



P(72 >T2\r2) = e- 



.i(2T/a) 



(27) 



With ( [27l i. the probability Ppsuc can be written as 
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Based on O and ([8]), and following a similar 

the probability 
where 



derivation as that for 



psuc 



can eventually be derived as 
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Putting (l23T l and (|29] | into ©, the coverage probability can 
then be directly obtained. 

V. Numerical results and discussion 

In this section, computer simulations are conducted to verify 
the correctness of the analytical results and investigate the 
performance of the two-tiers heterogenous cellular system. 
Simulation parameters are set as shown in Table I. 

Fig. 4 shows the coverage probability versus the power 
offset A. It is clear that the analytical results match well with 
the simulation results, demonstrating the correctness of the 
analytical results. Furthermore, the coverage probability in- 
creases with the power offset and achieves its maximum when 
the power offset equals to Pi/P-2, i.e., 20dB. The coverage 
probability drops slightly when the power offset is larger than 
20dB. It means that in order to achieve the maximum uplink 



TABLE I 
SIMULATION SETTINGS 



A2 


3.06 X 10"^ 


Ao 


7.66 X 10-^ 


Ti 


1 


T2 


1 


d 


50/\/3 


P1/P2 


100 



coverage probability, the power offset is suggested to be set 
equal to the ratio between the transmission powers of macro 
and pico BSs. This is not surprised since under this setting, 
the pico BSs deployed for coverage improvement can serve 
the mobile users with equal probability as the macro BS. 

The effect of the intensity of pico BSs on the coverage 
probability is next investigated and the results are shown in 
Fig. 5. Here the x axis in the figure denotes the intensity of 
pico BSs normalized by A2 shown in Table I. It is shown that 
the uplink coverage probability increases with the intensity of 
pico BSs as expected. This result demonstrates the benefit of 
pico deployment for the cellular coverage improvement. 




VI. Conclusions 

In this paper, an analytical approach to performance analysis 
for heterogenous cellular systems was proposed based on the 
theory of stochastic geometry. The uplink coverage probability 
was derived in closed-form and its correctness was verified by 
simulation results. From the analytical results, the impacts of 
both pico deployment and range expansion on the system's 
uplink performance were investigated. It has been found that 
in order to achieve the maximum uplink coverage probability, 
the power offset should be set equal to the ratio between the 
powers of macro and pico BSs. Meanwhile, the results have 
theoretically justified the benefit of the deployment of pico 
BSs. 
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